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We determine the variance-optimal hedge when the logarithm of 
the underlying price follows a process with stationary independent 



|Xy/| increments in discrete or continuous time. Although the general so- 

lution to this problem is known as backward recursion or backward 
• ■ stochastic differential equation, we show that for this class of pro- 

lan ' cesses the optimal endowment and strategy can be expressed more 

, explicitly. The corresponding formulas involve the moment, respec- 

tively, cumulant generating function of the underlying process and a 
Laplace- or Fourier-type representation of the contingent claim. An 
example illustrates that our formulas are fast and easy to evaluate 
numerically. 

> 

1. Introduction. A basic problem in mathematical finance is how an 
option writer can hedge her risk by trading only in the underlying asset. 
' This question is well understood in frictionless complete markets. It suffices 

\& . to buy the replicating portfolio in order to completely offset the risk. This 

elegant approach works well in the standard Black-Scholes or Cox-Ross- 
Rubinstein setup, but not much beyond. 
£3 . On the other hand, it has often been reported that real market data ex- 

hibit heavy tails and volatility clustering. Two common ways to account for 
such phenomena are some sort of stochastic volatility or jump processes, or a 
combination of both. In this paper we adopt the second approach and assume 
^ ' that the logarithmic stock price follows a general process with stationary, 

independent increments, either in discrete or continuous time. Processes of 
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this type play by now an important role in the modeling of financial data 
(cf. [3, 13, 14, 36]). 

Since replicating portfolios typically do not exist in such incomplete mar- 
kets, one has to choose alternative criteria for reasonable hedging strategies. 
If you want to be as safe as in the complete case, you should invest in a 
superhedging strategy (cf., e.g., [16]). In this case you may "suffer" profits 
but no losses at maturity of the derivative, which is very agreeable. On the 
other hand, even for simple European call options, only trivial superhedging 
strategies exist in a number of reasonable market models ( "buy the stock" ; 
cf. [8, 12, 22]). 

Alternatively, you may maximize some expected utility among all port- 
folios that differ only in the underlying asset and have a fixed position in 
the contingent claim. Variations of this approach have been investigated in 
[8, 9, 19, 30, 31]. 

In this paper we follow a third popular suggestion, namely, to minimize 
some form of quadratic risk (cf. [11, 21, 50, 52] for an overview). This can be 
interpreted as a special case of the second approach if we allow for quadratic 
utility functions. 

Quadratic hedging comes about in two different flavors: local risk min- 
imization as in [20] and [49], and global risk minimization (i.e., variance- 
optimal hedging, mean-variance hedging) as in [11] and [50]. Roughly speak- 
ing, one may say that locally optimal strategies are relatively easy to com- 
pute but hard to interpret economically, whereas the opposite is true for the 
globally optimal hedge. This paper focuses on the second problem, but as a 
by-product, we also obtain the locally optimal Follmer-Schweizer hedge. In 
discounted terms, the global problem can be stated as follows: If H denotes 
the payoff of the option and S denotes the underlying price process, try to 
minimize the squared L 2 -distance 

(1) E{{c + G T {0)-H) 2 ) 

over all initial endowments c£l and all in some sense admissible trading 
strategies 6. Here, G T (0) = Jq 9 t dS t [resp. G T {0) = £Li 9 n /\S n in discrete 
time] denotes the cumulative gains from trade up to time T. The idea is 
obviously to approximate the claim as closely as possible in an L 2 sense. 
Even though one may argue that one should not punish gains, the clarity 
and simplicity of this criterion are certainly appealing. Since it is harder to 
explain, we do not discuss local risk minimization here, but refer instead to 
[52]. 

By way of duality, quadratic optimization problems are related to (gener- 
ally signed) martingale measures, namely the Follmer-Schweizer or minimal 
martingale measure for local optimization, and the variance- optimal martin- 
gale measure for global optimization. A similar duality has been established 
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and exploited in many recent papers on related problems of utility maximiza- 
tion or portfolio optimization (cf. [7, 8, 17, 18, 23, 24, 25, 26, 33, 34, 40, 47]). 
Roughly speaking, the minimal martingale measure is the martingale mea- 
sure whose density can be written as 1 + Jq OtdMt for some 9, where M 
denotes the martingale part in the Doob-Meyer decomposition of S. The 
integrand 9 can be determined relatively easily in terms of the local behavior 
of S, which may be given by a stochastic differential equation or by one-step 
transition probabilities in discrete time. By contrast, the variance-optimal 
martingale measure is characterized by a density of the form c + J 9' t dSt 
for some c € R and some (generally different) integrand 9'. Here, it is usu- 
ally much harder to determine 9' . This holds with one notable exception, 
namely if the so-called mean-variance trade-off process is deterministic, in 
which case both measures coincide. More specifically, the integrands 9 and 
9' above tally because the difference Jq 9t dSt — Jq 9t dMt is a constant and 
can be moved to c. In this case of deterministic mean-variance trade-off, 
globally risk-minimizing hedging strategies can be computed from locally 
risk-minimizing ones. The setup in this paper is among the few models of 
practical importance where the condition of deterministic mean-variance 
trade-off naturally holds. 

The process formed by conditional expectation of the option's payoff un- 
der the minimal, respectively, variance-optimal martingale measure is some- 
times interpreted as a derivative price process. In jump-type models one 
has to be careful at this point because these measures are generally signed 
and may lead to arbitrage. In particular, the variance-optimal initial capital 
for a positive payoff may be negative; see [2], Section 4 and the numerical 
example at the end of Section 5.2 below. 

Even in the case of deterministic mean- variance trade-off, the actual com- 
putation of variance-optimal hedging strategies involves the joint predictable 
covariation of the option's "price process" and the underlying stock. For 
general claims, how to obtain this covariation may not seem evident. It can 
be computed quite easily if the payoff is of exponential type e zXr , where 
X := log(5/5o) denotes the process with stationary, independent increments 
driving the stock price S. The reason is that the "price process" for such 
exponential payoffs under the variance-optimal martingale measure is again 
the exponential of a process with stationary independent increments, which 
leads to handy formulas for the corresponding hedge. Since the optimality 
criterion in (1) is based on an L 2 -distance, the resulting hedging strategy 
is linear in the option. This suggests writing an arbitrary claim as a lin- 
ear combination of exponential payoffs. Put differently, we work with the 
inverse Laplace (or Fourier) transform of the option. This will be done in 
Section 2 for discrete-time and in Section 3 for continuous-time processes, 
respectively. One could go even one step further and generalize the results 
to arbitrary processes with independent increments because they still share 
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the important property of deterministic mean-variance trade-off. However, 
we chose not to do so in order not to drown the arguments in technicalities 
and because this more general class plays a minor role in applications. 

Since the first version of this paper circulated, the idea of using Fourier 
or Laplace transforms with Levy processes has been applied independently 
in the framework of option pricing by [4] as well as [42] and, very recently, 
in the context of quadratic hedging by [53] . Explicit integral representations 
of a number of concrete payoffs are to be found in Section 4. 

Section 5 illustrates the application of the results. We compare the vari- 
ance-optimal hedge of a European call in a pure-jump Levy process model to 
the Black-Scholes hedge as a benchmark. Since the results in the subsequent 
sections rely heavily on bilateral Laplace transforms, the Appendix contains 
a short summary of important results in this context. 

To keep the presentation and notation simple, we confine ourselves to one 
single underlying asset. Extensions to the multivariate case and to path- 
dependent claims will be provided elsewhere. For unexplained notation, we 
refer the reader to standard textbooks on stochastic calculus such as, for 
example, [29] or [41]. 

2. Discrete time. Let (fl, J 7 , (^Fn)ne{o,i,...,N}-> P) denote a filtered prob- 
ability space and let X = (Jf n )n=o,i,...,7V denote a real-valued process with 
stationary, independent increments in the sense that: 

1. X is adapted to the filtration (J~ n )ne{0 l N}> 

2. X = 0; 

3. AX n := X n — X n _i has the same distribution for n = 1, . . . ,N; 

4. AX n is independent of T n -\ for n = 1, . . . , N. 

We consider a nondividend-paying stock whose discounted price process S 
is of the form 



with some constant So > 0. We assume that E(Si) = S$E(e 2Xl ) < oo, which 
implies that the moment generating function m:z^ E(e zXl ) is defined at 
least for z € C with < 5R(z) < 2. Moreover, we exclude the degenerate case 
that S is deterministic. Put differently, Var(e^ 1 ) = m(2) — m(l) 2 does not 
vanish, which can be viewed as a no arbitrage condition. 

Our goal is to determine the variance-optimal hedge for a European- 
style contingent claim on the stock that expires at N with discounted pay- 
off H. Mathematically, H denotes a square-integrable, .F/v-nieasurable ran- 
dom variable of the form H = /(Sn) for some function / : (0, oo) — ► R. More 
specifically, we assume that / is of the form 



S n = Sbexp(X n ) 



(2) 
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for some finite complex measure II on a strip {z G C : R' < 5ft(z) < R}, where 
R',ReR are chosen such that E(e 2R ' Xl ) < oo and E(e 2RXl ) < oo. The 
integral representation of European call and put options, and several other 
payoffs can be found in Section 4. 

Remark 2.1. Loosely speaking, the option's payoff at maturity is writ- 
ten as a linear combination of powers of the S or exponentials of X. Put 
differently, its payoff function is a kind of inverse Mellin or Laplace trans- 
form of the measure LI. To be more specific, let us consider the case R' = R, 
that is, II is concentrated on the straight line R + iWL. Denote by £ the inverse 
Laplace transform of IL in the sense that £(x) = f^i™ e zx U(dz) for x G R. 
Then 



Up to a factor e , the function £ is just the Fourier transform of a complex 
measure on the real line [namely the measure v with v{B) = Ii{R + iB) for 
Borel sets BcK]. The reason to allow also R ^ is simply that £ cannot be 
written as the Fourier transform of a finite measure for important claims as, 
for example, European calls. Examples with R' < R are given in Remark 4.1 
part 2. 

The variance-optimal hedge minimizes the L 2 -distance between the op- 
tion's payoff and the terminal value of the hedging portfolio. To be more 
specific, define the set of admissible strategies as the set of all predictable 
processes 9 such that the cumulative gains G n (9) := YX=i Ok^Sk are square- 
integrable for n= 1, . . . , N. We call cf) G a variance- optimal hedging strategy 
and Vq G R a variance- optimal initial capital if c = Vo and 6 = <p minimize 
the expected squared hedging error 



over all initial endowments cGK and all admissible strategies 9 G ©. Let 
us emphasize that the variance-optimal initial capital is, in general, not an 
arbitrage-free price and can be negative for a positive payoff; see Section 5.2 
for a concrete example. 

In our framework the variance-optimal hedge and its corresponding hedg- 
ing error can be determined quite explicitly: 

Theorem 2.1. The variance-optimal initial capital Vq and the variance- 
optimal hedging strategy 4> ar ^ given by 



H = f(S N ) = /(exp(A0v + log(So))) = i(X N + log(5 )). 



(3) 



E((c + G N (9)-H) 2 ) 




V = H 



(5) 



(f>n = in + 



A 



{Hn-x-Vo-Gn-xiffy) 
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where the processes H n and £ n and the constant A are defined by 
. m(z + 1) — m(l)m(z) 
m{2) — m(l) 2 

h(z) := m(z) - (m(l) - l)g(z), 
(6) A — 



m(2)-2m(l) + l' 

T/ie optimal hedge Vq,4> is unique up to a null set. 

Remark 2.2. One may also consider a similar problem where the initial 
endowment is fixed and the mean squared difference in (3) is minimized only 
over the strategies 9 E 0. This risk-minimizing hedging strategy for given 
initial capital c is determined as in Theorem 2.1 with Vq = c instead of (4). 

Theorem 2.2. The variance of the hedging error E({Vq + Gn{4>) ~ H) 2 ) 
in Theorem 2.1 equals 



Jo-= J J Jo( yi z)U(dy)U(dz), 
m(2)-m(l) 2 



where 

a(y,z) ■■=Hy)Hz, m{2) _ 2mil) + 1 , 

b(y,z) :=m(y + z) 

— (m(2)m(y)m(z) — m(l)m(y + l)m(z) 

— m{l)m(y)m{z + 1) + m(y + l)m(z + l))(m(2) — m(l) 2 )" 1 , 

{ c y+z u( My,z) N -m(y + z) N . . 

So b(y,z) r , ifa{y,z)^m{y + z), 

a(y,z) -m(y + z) 

So +z b(y, z)Nm{y + z) N ^, if a(y, z) = m{y + z). 

The remainder of this subsection is devoted to the proofs of Theorems 
2.1 and 2.2. As noted by [51], the variance-optimal hedge can be obtained 
from the option's Follmer-Schweizer decomposition if the so-called mean- 
variance trade-off process of the option is deterministic. The latter is defined 
as 

" (EjAS^.,)) 2 (m(l)-l) 2 
n - ^ Var(A5 fc [^ fc _i) ?n(2)-m(l) 2 ' 



VARIANCE-OPTIMAL HEDGING 



7 



The Follmer-Schweizer decomposition plays a key role in the determination 
of locally risk-minimizing strategies in the sense of [20] and [49], and it is 
defined as follows. 

Definition 2.1. Denote by S = So + M + A the Doob decomposi- 
tion of S into a martingale M and a predictable process A. The sum 
H = Hq + J2n=i £,nAS n + Ln is called a Follmer-Schweizer decomposition 
of H € L 2 (P) if is .^-measurable, £ € © and L is a square-integrable 
martingale with Lq = that is orthogonal to M (in the sense that LM is a 
martingale). We will use this notion as well if H,Hq,^ and L are complex- 
valued, in which case we require € and £ 0. 

In discrete time any square-integrable random variable admits such a de- 
composition, which can be found by a backward recursion (cf. [51], Proposi- 
tion 2.6). However, since this method does not yield a closed-form solution 
in our framework, we do not use these results. Instead we proceed in two 
steps. First, we determine the Follmer-Schweizer decomposition for options 
whose payoff is of power type. Second, we consider claims which are lin- 
ear combinations of such options in the sense of (2). Here, we rely on the 
linearity of the Follmer-Schweizer decomposition in the claim. 

Lemma 2.1. Let z e C with Sf e L 2 (P). Then H(z) = S Z N admits a 
Follmer-Schweizer decomposition H(z) = H {z)q + ^2^ =1 ^{z) n /S.S n + L{z) N , 
where 

H(z) n = h(z) N ~ n S z n , 

(7) i{z) n = g{z)h{z) N - n S z n Zl 

n 

L(z) n = H{z) n - H(z) - £ C(z) k AS k 

k=l 

and g(z),h(z) are defined in Theorem 2.1. 

Proof. The statement could be derived from Proposition 2.6 and Lemma 
2.7 of [51], but it is easier to prove it directly. Since Sf is square-integrable, 
all the involved expressions are well defined. From (7) it follows that 

(8) AL(z) n = S^Mz) N - n (e zAX " - h{z) - g(z)(e AX " - 1)). 
Since 

E{e zAXn - h{z) - g(z)(e AXn - 1)) 

(9) = m{z) - h{z) - g{z){m{l) - I) 

= 0, 
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this implies that E(AL(z) n \J- n -i) = and hence L(z) is a martingale. 
The Doob decomposition S = So + M + AoiS satisfies 

(10) AA n = E{AS n \T n ^) = S n _i(m(l) - 1) 

and hence AM n = S n ^\(e AXn — m(l)). In view of (8), we obtain 

AM n AL(z) n = S*±\h(z) N - n (e AXn -m(l))(e zAX " -h{z) - g{z){e AXn - 1)). 

From 

E( e AXn (e zAXn - h(z) - g(z)(e AXn - 1))) 

= m(z + 1) — h(z)m(l) — g(z)m(2) + g(z)m(\) 
= 

and (9) it follows that E(AM n AL(z) n \J r n -i) = and hence ML(z) is a 
martingale as well. □ 

Proposition 2.1. Any contingent claim H = /(SV) as in the begin- 
ning of this subsection admits a Follmer-Schweizer decomposition H = Hq + 
2~2n=i £nAS n + Ljv- Using the notation of the previous lemma, it is given by 

H n = J H(z) n U(dz), 

Cn = Jt(z) n U(dz), 

/n 
L(z) n U(dz) =H n -H -J2 tk&S k . 
k=i 

Moreover, the processes H n , £ n and L n are real-valued. 

PROOF. First, note that / E(\AL(z) n \ 2 )\U\(dz) < oo, where | TT | denotes 
the total variation measure of II in the sense of [44], Section 6.1. From 
Fubini's theorem we conclude that 

E(AL n l B ) = E^J AL(z) n U(dz)l B ^j 

= J E(AL(z) n l B )U(dz) = 

for any B £ J~ n -i- Hence L is a martingale. Similarly, it is shown that ML 
is a martingale as well. The assertion concerning the decomposition follows 
from Lemma 2.1. 

Since H and S n are real-valued, we have 

TV 

= (Hq - Hq) + UAS n + (L N - L N ), 

ii=i 
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which implies that = 3f(# ) +E»=i 3 ! (?n)A5 n + 3(L 7V ). Since the Follmer- 
Schweizer decomposition of is unique (cf. [48], Remark 4.11), we have that 
Hq, £ n and L n are real- valued for n = 1, . . . , N. □ 

Finally, we apply the preceding results to determine the variance-optimal 
hedge. 

Proof of Theorem 2.1. As observed by [48], Proposition 5.5, the 
process S has deterministic mean- variance trade-off. From Proposition 2.1 
and [51], Theorem 4.4, it follows that the variance-optimal hedging strategy 
<f> satisfies 

<t>n=£,n + K(H n -l - Hq- G n _l(0)) 

with 

AA n X 

[cf. (10)]. Moreover, the variance-optimal initial capital equals Vo- 

For the uniqueness statement, suppose that Vo G JR and <f> G G lead to a 
variance-optimal hedge as well. Define Vo := \{Vq + Vo) and (f> := ^((f> + 4>) G 
Q. It is easy to verify that we would have 

E((V + G N (4>) - Hf) < lE((V + G N (4>) - H) 2 ) 

+ \E{(y Q + G N Q>)-H) 2 ), 

contradicting the optimality of (Vo,<^) and (Vo,</>), if Vo + Gjy(4>) and Vo + 
Gn{4>) did not coincide almost surely. Hence 

V o -Vo + G N (4>-0) = O. 

In particular, Gn{<P — <f) is J~n-\ -measurable. We obtain 

O = Var(GW(</>-0)|^v-i) 

= Var((0 - 4>) N /S.S N \T N -i) 

= ((0-0) N S N - l ) 2 (m(2)-m(l) 2 ), 

which implies that (cp — 4>)n = almost surely. By induction, we conclude 
that (<f) — 4>) n = for n = N — 1, . . . , 1 and hence also Vq = Vq. 

The remark following Theorem 2.1 follows from [51], Proposition 4.3. □ 

Proof of Theorem 2.2. According to [51], Theorem 4.4, the variance 
of the hedging error equals 

N N 

(ii) ]T£((AL n ) 2 ) n a-v^) 

n=l k=n+l 
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with Afc = \/Sk-i and AA k as in (10). Since AL n = J AL(z) n TL(dz), we have 
that 

(AL„) 2 = J J AL(y) n AL(z) n U(dy)n(dz) 

and hence 

(12) E({AL n f) = J J E(AL(y) n AL(z) n )U(dy)U(dz) 

by Fubini's theorem. Equation (8) implies 
AL(y) n AL(z) n = S^hM"^)"^**** - h(y) - g(y)(e AX " - 1)) 

x (e zAX ™ - h(z) - g(z)(e AXn - 1)). 
bince E(Sltl) = Sl +Z m(y + z)"" 1 and so forth, we have 

E(AL(y) n AL{z) n ) = S y +Z (h{y)h{z)) N ~ n m(y + z) n - l b(y, z) 

with 

b(y, z) = m{y + z) - m(y)(h(z) - g{z)) - m(y + l)g(z) 

- (h(y) - g(y))(m(z) - h(z) + g(z) - g(z)m(l)) 

- g(y)(m(z + 1) - m(l)(h(z) - g(z)) - g(z)m(2)). 

This expression coincides actually with b(y, z) in the statement of the theo- 
rem. Consequently, we have shown 

N N 

-£E(AL(y) n AL(z) n ) ]J (l-X k AA k ) 

n=l k=n+l 

= S y+ Z b(y, z)a(y, z) N ^ £ (^f^) ^ 

f^ x \ a(y,z) J 

- ^y +z h(„ Mv, z ) N ~ m (y + z) N 

a{y,z) -m(y + z) 

unless the denominator vanishes in the last equation. In view of (11) and 
(12), we are done. □ 

Let us briefly discuss the structure of the variance-optimal hedge. The 
process £ in the Follmer-Schweizer decomposition coincides with the locally 
risk-minimizing strategy. The process H n = Hq + Ylk=i Zk^Sk + L n that 
appears on the right-hand side of the Follmer-Schweizer decomposition may 
be interpreted as a "price process" of the option. However, since this process 
may generate arbitrage, one should be careful with this interpretation. Note, 
however, that the difference between the locally and globally optimal hedging 
strategy in (5) is proportionate to the difference between this "option price" 
H n —\ and the investor's current wealth. 
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3. Continuous time. We turn now to the continuous-time counterpart 
of the previous section. Similarly as before, (Q, J 7 , (Ft)te[o,T\ ; P) denotes a 
filtered probability space and X = {Xt) te r 0T i denotes a real- valued process 
with stationary, independent increments (PUS, Levy process) in the sense 
that: 

1. X is adapted to the filtration (J~t)t€[oT\ an d has cadlag paths; 

2. X = 0; 

3. the distribution of X t — X u depends only on t — u for < u < t < T; 

4. Xt — X u is independent of T u for < u < t < T. 

As in the discrete-time case, the distribution of the whole process X is de- 
termined by the law of X±. The latter is an infinitely divisible distribution 
which can be expressed in terms of its Levy-Khinchine representation. Al- 
ternatively, one may characterize it by its cumulant generating function, 
that is, by the continuous mapping k:D^C with E(e zXt ) = e tK ^ for 
z € D := {z € C : E(e ft ^ Xl ) < oo} and t £ M + . For details on Levy processes 
and unexplained notation, we refer the reader to [29, 41, 46]. 

The discounted price process S of the nondividend-paying stock under 
consideration is supposed to be of the form 

S t = S exp(X t ) 

with some constant 5o > 0. Again, we assume that E(Sf) = SQE(e 2Xl ) < 
co, which means that z & D for any complex number z with < 3t(z) < 2. 
Moreover, we exclude the degenerate case that S is deterministic, that is, 
we have k(2) — 2k(1) ^ 0. This can be viewed as a no arbitrage condition. 

As in Section 2 we consider an option with discounted payoff H = /(St), 
where / is given in terms of a finite complex measure II [cf. (2)]. The choice 
of the set of admissible trading strategies is a delicate point in continuous 
time. Following [50], Section 1, we choose the set 

8:= |^€L(5):^ 9 t dS t GH 2 |, 

which is well suited for quadratic optimization problems. Here, the space 
TL 2 of semimartingales is defined as follows: 

Definition 3.1. For any real- valued special semimartingale Y with 
canonical decomposition Y = Yq + N + B, we define 

[|Y|| W 2 := [[Yolb + \\y/[N, N] T \\ 2 + || var(B) T || 2 , 

where var(£>) denotes the variation process of B and || • [[2 denotes the 
L 2 -norm. By 7i 2 we denote the set of all real- valued special semimartingales 
Y with ||y||-^2 < 00. 
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In our setup, this set can be expressed more easily as follows: 
Lemma 3.1. We have 

6 = i.9 predictable process: ^(^J ^t^t- dtj < oo|. 

Proof. From Lemma 3.2 below we conclude that At = k(1) J S u -du 
and 

(13) (M,M) t = (K(2)-2«(l)) f Sl_du 

Jo 

for the canonical decomposition 5 = Sq + M + A of the special semimartin- 
gale S. Hence we have 

(14) A t = f X u d(M,M) u 

Jo 

with A u := \/S u ~ and A := k ( 2 )^2k(i) 1 Therefore, the mean- variance trade-off 
process 

in the sense of [50], Section 1, is deterministic and bounded. According 
to [50], Lemma 2, we have that 9 £ holds if and only if 9 is predictable 
and E(fu9?d{M,M)t) < 00. Since 

■ d(M, M) t = («(2) - 2*c(l)) l T 9 2 t S^ dt, 



t 

Jo 
the assertion follows. □ 

If we denote by Gt(9) := ^9 u dS u the cumulative gains process of 9 € 
0, then the variance- optimal initial capital and variance- optimal hedging 
strategy can be defined as in the previous section (with T instead of N). 

The following characterizations of the variance-optimal hedge and its ex- 
pected squared error correspond to Theorems 2.1 and 2.2. They are proved 
at the end of this subsection. 

Theorem 3.1. The variance-optimal initial capital Vq and the variance- 
optimal hedging strategy 4> are given by 

V = H 

and the expression 

(15) <h = & + 4~( H *- " y o " G *-(^))> 
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where the processes Ht and £t and the constant A are defined by 

_ k(z + 1) — k(z) — k(1) 
1[Z) : ~ k(2)-2k(1) ' 

rj(z) := k(z) - k(1)j(z), 
«(1) 



(16) A 



k(2) - 2«(1)' 
fli : = J Sfe^ T -^Il(dz), 



JTie optimal initial capital is unique. The optimal hedging strategy (j>t(w) is 
unique up to some (P(duj) <8> (it) null set. 

Remark 2.2 on risk-minimizing hedging for fixed initial endowment c ap- 
plies in continuous time as well. 

Theorem 3.2. The variance of the hedging error E{{Vq + Gt{4>) — H) 2 ) 
in Theorem 3.1 equals 



Jo--=f J J (y,z)u(dy)n{dz), 



where 



a(y,z) :=rj(y) +r/(z 



k(2)-2ac(1)' 

P(y,z) := K(y + z) - n(y) - k{z) 

( K (y + 1) - K(y) - k(1))(k(z + 1) - k{z) - «(!)) 
k(2)-2k(1) 

{ e a(y,z)T _ e K{y+z)T 

St z p(y,z)— r — ■ — -, ifa(y,z)^K(y + z), 
a{y,z) - K{y + z) 

S y +Z p(y, z)Te K ^ T , if a(y, z) = n{y + z) . 

Remark 3.1. If (fj,, cr 2 , v) denotes the Levy-Khinchine triplet of X (rel- 
ative to the truncation function x i— > ^l-ru^n), then we have 

a 2 f 

k(z) =fiz + —z 2 + / {e zx - 1 - zxt { \ x \< 1} )v(dx) 

for z G D (cf. [46], Theorem 25.17). In particular, we have k(z) = \iz + 
(a 2 /2)z 2 for Brownian motion. Note that 



I rR+ioo e (x-fi)z+(a 2 /2)z 2 



=— / Cb 
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for any R > 0, where $ denotes the cumulative distribution function of 
N(0, 1). Using the same decomposition and substitution as in Remark 4.1, 
one easily shows that Vq and (ft in Theorem 3.1 coincide with the Black- 
Scholes price and the replicating strategy in the case of a European call H 
and Brownian motion X. This does not come as a surprise because perfect 
hedging is clearly variance-optimal. 

The hedging strategy in Theorem 3.1 is given in feedback form, that 
is, it is only known in terms of its own gains from trade up to time t. 
From a practical point of view, these gains are obviously known to the 
trader. However, they cannot be computed recursively as in the discrete- 
time case. Therefore, one may prefer an explicit expression for Gt(4>) from 
a mathematical point of view. It is provided by the following theorem. 

Theorem 3.3. Suppose that P(AX t = log(l + l/A) for some te (0,T]) = 
0. Then the gains process of the variance-optimal hedging strategy <f> in The- 
orem 3.1 is of the form 



Remark 3.2. The condition on X is equivalent to assuming that the 
Levy measure of X puts no mass on log(l + 1/A). This holds for any model 
of practical importance. In the general case it is still possible to give an 
explicit, albeit more involved, representation of the gains process (cf. [28], 



Moreover, observe that X and Y are both Levy processes (cf. [32], Lemma 
2.7 and straightforward calculations). Recall that the stochastic exponential 
£(U) of a real-valued Levy process or any other semimartingale U can be 
written explicitly as 




where the processes X and Y are defined as 



(17) 





(6.8)). 



£(U) t = exp(U t - \[U, U) t ) J] (1 + MJ U ) exp(-AC4 + ±(A[/ U ) 2 ) 



u<t 



(cf. [41], Theorem 11.36). 
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The remainder of this section is devoted to the proof of Theorems 3.1-3.3. 
The approach parallels the one in the previous section. As before, we deter- 
mine the Follmer-Schweizer decomposition of the claim and apply results 
that relate this decomposition to the variance-optimal hedge. 

Lemma 3.2. Let z eC with Sf, e L 2 (P). Then S z is a special semi- 
martingale whose canonical decomposition Sf = Sq + M(z) t + A(z) t satisfies 

(18) A{z). = k{z) C S z u _du 

Jo 

and 

(19) (M{z),M) t = ( K (z + l)- K {z)-K(l)) f S^du, 

Jo 

where M = M(l) corresponds to z = 1 as in the proof of Lemma 3.1. 

Proof. Note that almost by definition of the cumulant generating func- 
tion, N(z) t := e~ K ^ t Sl is a martingale. Integration by parts yields Sf = 
e K ^ t N(z) t = S^ + M(z) t + A{z\ with M(z) t = /* e R ^ s dN{z) u and A(z) as 
claimed. Moreover, we have 

[M(z),M] t = [S z ,S] t 

= S t ~^~ — Sq~^~ — / S u _dS u — / S u — dS u 
Jo Jo 

= M(z + 1), - f S z u _ dM u - f S u „ dM(z) u 
Jo Jo 

+ (k(z + 1) - k(z) - k(1)) /* Sit 1 du. 

Jo 

Note that the first three terms on the right-hand side are local martingales. 
Since (M(z),M) is the predictable part of finite variation of the special 
semimartingale [M(z),M], equation (19) follows. □ 

Definition 3.2. Denote by S = So + M + A the canonical special semi- 
martingale decomposition of 5 into a local martingale M and a predictable 
process of finite variation A. The sum H = Hq + Jq ^tdSt + Lt is called 
a Follmer-Schweizer decomposition of H € L 2 (P) if Hq is J-o-measurable, 
£ £ O and L is a square-integrable martingale with Lq = that is orthogonal 
to M (in the sense that LM is a local martingale) . We will use this notion as 
well if H, Hq, £ and L are complex- valued, in which case we require 5?(£) € 
and € 6. 
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The existence of a Follmer-Schweizer decomposition was established in [50], 
Theorem 15, in our case of bounded mean-variance trade-off. It can be ex- 
pressed in terms of a backward stochastic differential equation. Since the 
latter may be hard to solve, we do not use this result. Instead, we prove 
directly that the continuous-time limit of the expressions in Section 2 leads 
to a Follmer-Schweizer decomposition. 

Lemma 3.3. Let z E C with G L 2 (P). Then H(z) = Sj, admits a 
Follmer-Schweizer decomposition H{z) = H(z) + J T ^{z) t dS t + L(z) T , where 

H(z) t :=e^ T -^Sl, 

(20) ^(z) t :=l(z)e^ T ^S^\ 

L(z) t :=H(z) t -H(z) Q - f\{z) u dS u 

Jo 

and "f(z),r}(z) are defined in Theorem 3.1. Moreover, M is a square-integrable 
martingale and hence L(z)M is a martingale. 

Proof. Partial integration and (18) yield 

H(z) t = H(z) Q + C e^ { - T - u UM(z) u + {n{z) - n{z)) C e^ T ~ u ^S^ du 
Jo Jo 

and 

[\(z) u dS u = [\(z) u dM u + K(l) 7 (z) f e^ T ~^S z u _du. 
Jo Jo Jo 

Since k{z) — n(z) — k(1)j(z) = 0, the predictable part of finite variation in 
the special semimartingale decomposition of L(z) vanishes and we have 

(21) L(z) t = f e^ T - u UM(z) u - l\(z) u dM u , 

Jo Jo 

which implies that L{z) is a local martingale. 

From (19) for z and 1 instead of z it follows that 

(L(z\M) t = f 'e^ T ' u U(M(z),M) u - f\{z) u d{M,M) u 
Jo Jo 

= (k(z + 1) - k(z) - k(1) - 7(z)(k(2) - 2k(1))) 

x (* e^){T-u) s z+i du 
Jo 

= 0. 

Consequently, L(z)M is a local martingale as well. 
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Similar calculations yield 



(L(z),L(z)) t = (L(z),L(z)) t 

\k{z + 1) - k{z) - k(1)| 2 



(22) = (/c(23i(z)) - 2®.{k{z)) 



and 



(23) 



Since 



[ T \az) t \ 2 sidt- 

Jo 



k(2)-2k(1) 



k(z+ 1) - k(z) - k(1) 



k(2)-2k(1) 







£(S f 2R(z) ) = S(5f (z) ) = <?f (2) e^ 2S ^» < 5 2SR W(1 V e T < m ^) < oo, 

(24) 

it follows that E((L(z), L(z))t) < oo. Therefore, L is a square-integrable 
martingale. 

Similarly, (23) and (24) yield that 9?(f (2)) G 9 and 3 (£(*)) G 6. Equa- 
tions (19) and (24) for 1 instead of z imply that M is a square-integrable 
martingale. □ 

Lemma 3.4. There exist constants c\, . . . , C5 > such that 

(25) »fo(z)) < ci, 

o,^))- 2RW2 ))-Wi±|=^(* 

(26) <-c 2 5?(r/(z))+C3, 

|7(^)! 2 <-C43?(77(z))+ C5 

/or any z G C with R' < 3?(z) < i?, where 7 and 77 are defined as in Theorem 
3.1. 

Proof. Since «; is continuous, there is a constant c% > such that 

(27) l«(2»(z))| <2c 6 



for any z with i?' < K(z) < R. Since (L(z), L(z)) is increasing, (22) yields 

4 « w) - TO - Mi±|=gg_^ >o. 
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In particular 

»(k(z)) < |«(2»(z)) < c 6 

and 

/ \ \k(z + 1) — k(z) — k(1)\ 2 ml , 

(28) K (2)-L(l) S*»-«t«M). 

which implies 



| K (l) 7 (z)| 2 < C7 -c 8 K( K (z)) 
1 

4 



< C 2 + I(K( K (z))) 2 



<(|^( K (z))|+ C9 ) 



2 



for some C7, cs > and eg := y C7 + 4c§. This yields 

R(77(^)) = »(«(«)) -K(«(1) 7 W) 
<»(«(*)) + Klfr(*)| 

(29) 

< Ci + 

< eg + 2c6 =: ci 
with cio := eg + §C6- Inequality (28) also yields 

h(z)\ 2 < Cll - C ^(K(z)) 

for some cn,C4 > 0, which, together with (29), leads to 

| 7 (^)| 2 < en - c 4 (K(r ? (^)) - cio) = c 5 - c 4 5?(r7(z)) 

with C5 := en +C4C10. Finally, the second inequality in (26) follows from 
(27), (29) and k(2) - 2k(1) > 0. □ 

Proposition 3.1. Any contingent claim H = /(St) as in the begin- 
ning 0/ this subsection admits a Follmer-Schweizer decomposition H = Hq + 
Jq £t dSt + Lt ■ Using the notation 0/ Lemma 3.3, it is given by 

H t = J H(z) t U(dz), 

(30) e* = / £(*) t n(<te), 

it = y L(z) t n(cfe) = H t -H -J* i u dS u . 
Moreover, the processes Ht, £t an d Lt are real-valued. 
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Proof. Let z e C with R' < 3t(z) < Since \H(z) t \ 2 = ^"WK 1 "-*) sf* w , 



we have that (-z)t | 2 ) is bounded by some constant which depends only 
on R and R' [cf. (24) and (25)]. It follows that H t is a well-defined square- 
integrable random variable. Similarly, (22), (24) and Lemma 3.4 yield, after 
straightforward calculations, that 



is bounded as well by such a constant. Therefore, Lt is a well-defined square- 
integrable random variable as well. Finally, (23) and Lemma 3.4 yield that 
E(\i(z) t S t -\ 2 ) and also E(J^ \£{z) u \ 2 Sl_ du) are bounded by some constant 
which depends only on t, R and R' . Therefore, £ is well defined, and 3ft (£) € @ 
and € O by Lemma 3.1. The same Fubini-type argument as in discrete 
time shows that E((L t -L u )l B ) =0 and E((M t L t - M u L u )t B ) = for u < t 
and B S T u (cf. Proposition 2.1). Hence, L is a square-integrable martingale 
which is orthogonal to M. To be precise, we interpret L as the up to in- 
distinguishability unique modification whose paths are almost surely cadlag 
(cf. [41], Corollary LI). By Fubini's theorem for stochastic integrals (cf. [41], 
Theorem IV. 46), we have 



Together with (30) and (20) it follows that Hq,^ and L do indeed provide 
a Follmer-Schweizer decomposition of H . As in the proof of Proposition 
2.1, this time using [38], Theorem 3.4, instead of [48], the uniqueness of the 
real- valued Follmer-Schweizer decomposition yields that the processes Ht, 
£i and L t are real-valued. □ 

PROOF of Theorem 3.1. According to the proof of Lemma 3.1, the 
mean- variance trade-off process of S in the sense of [51], Section 1, equals 



In view of Proposition 3.1, the optimality follows from Theorem 3 and Corol- 
lary 10 of [50]. 

As in the proof of Theorem 2.1, it follows that Vo = Vo and Gt(<P) = 
Gt(<P) if Vo and (j> denote another variance-optimal hedge. Observe that the 
local martingale N t := - ^\ u dM u satisfies {N,N) T = ^ ' \ 2 u d(M,M) u = 
Kt, where X u is defined as in the proof of Lemma 3.1. From [5], Propositions 
3.7 and 3.9(h), and the remark after Definition 5.4, it follows that G{4> — <f)) 
is a £(N) martingale in the sense of that paper. By Proposition 3.12(i) in 



E(\L(z) t \ 2 ) = E({L(z),L(z)) t )<E({L(z),L(z)) T ) 
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the same paper, it is determined by its terminal value Gt{4> — = 0, that 
is, Gt{4> — (ft) = for any t € [0, T]. Hence 

O = £([G(0-0),G(0-£)] T ) 

i-T 



E^j\^-4>) 2 t d[M,M])j 
E^j\^-^) 2 t d(M,M)^ 
(k(2) - 2k(1))E^ (eft - 4>f t Sl dt^j . 



This implies that <fo(u;) = 4>t(u) outside some (P(doj) <X> <ii) null set. □ 

Proof of Theorem 3.2. Similarly as in Lemma 3.2, it is shown that 

(M(y),M(z)) t = ( K (y + z) - K (y) - k(z)) f S y u t z du. 

Jo 

From (21), (L(y),M) = and (19), it follows that 

(L(y),L(z)) t = f e (v(y)+v(z))(T-u) d{M{y):M{z))u 
Jo 

(31) - / < 7 (z) e (^)+^))( T -«) S z u Z l d(M(y),M) u 

Jo 

= (3(y, z) f e ^y^(z))(T-u) gy+z ^ 
Jo 

Consequently, 

(32) [ T e-VT-KJ d(L(y),L(z)) t = 0(y, z) S^ z e a ^ T ^ dt, 



n 



where K denotes the mean-variance trade-off process as in the proof of 
Lemma 3.1. Since E(S^ Z ) = SQ +z e K ( y+z ^ t , an application of Fubini's theo- 
rem yields 



E' 



f e -(KT-Kt) d ( L ( y ^ L ( z ^ =SV +Z /3(y,z) ^ e <V+z)t+a{y,z){T-t) ^ 



which equals Jo(y,z) 
Observe that 



ft(L(y),L(z)) = U(L(y) + m,L(y) + L(z)) - (L(y),L(y)) - (L(z),L(z))) 
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and 

(L(y) + L(z),L(y) + L(z)) 

< (L(y) + L(z),L(y)+L(z)) + (L(y) - L(z), L(y) - L(z)) 

= 2(L(y),L^)) + 2(L(z),L(z)). 

In the proof of Proposition 3.1 we noted that E({L(z) , L(z))t) and hence 
also the expected total variation of 5?((L(y), L(z))t) is bounded by some 
constant which depends only on R and R' . By replacing L(z) with iL(~z), it 
follows analogously that the total variation of 9((L(y), L(z))t) is bounded 
by a similar constant. Therefore, 

J J(L(y),L(z)) t n(dy)U(dz) 

is a well-defined continuous, predictable, complex-valued process of finite 
variation. Since 

L 2 t =J J L(y) t L(z) t U(dy)U(dz), 

an application of Fubini's theorem yields that 

L\- J J(L(y),L(z)) t U(dy)n(dz) 

is a martingale. This implies 

(L,L) t = J J(L(y),L(z)) t U(dy)n(dz) 

by definition of the predictable quadratic variation. Another application of 
Fubini's theorem yields 

j\-(K T -K t ) d{LjL)t = J J £ e -(KT-K t ) d{L{y ^ L{z))tU{dy)u{dz) 

and hence 

E {[ e ~ [KT ~ Kt)d ^ L )) 

= J J E^£e^ K ^d(L(y),L(z)) t y(dy)U(dz) 

= JI My,z)U(dy)U(dz). 

By [50], Corollary 9, the left-hand side of the previous equation equals the 
variance of the hedging error. □ 

Finally, we prove the explicit representation of the gains process. 
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PROOF of Theorem 3.3. By (15), G(<p) solves the stochastic differen- 
tial equation 

GM= t(^+ XiHu -~ Vo) )dS u - f^G u ^)dS u 



V Ju- J JO Ju- 

= [ (t; u S u - + \(H u --Vo))dX u + [ G s _(<f>)d(-\X) U . 
Jo Jo 

By [28], (6.8), this equation has a unique solution, which is given by 

Jo £{-\X) u 
Since £{— XX) U = (1 — \AX U )£(— AX) M _, the assertion follows. □ 

4. Integral representation of payoff functions. In the previous sections 
we used an integral representation of the payoff of the form (2). A precise 
characterization of the class of functions that allow for such a representation 
is delicate (cf. [6]) and of minor interest in this context. To derive hedging 
strategies and so forth, explicit formulas for II are required. They are pro- 
vided here for a number of payoffs. 

The basic example is, of course, the European call option H = (£V — K) + . 
Its integral representation (2) is provided by the following lemma. 

Lemma 4.1. Let K > 0. For arbitrary R > 1, s > 0, we have 

( s -K) + = — s z -^ -dz. 



2iri jR-ioo z(z - 1) 
PROOF. For K(z) > 1, we have 

/oo 
(e x -K) + e~ zx dx 
-oo 



K 



1-2 



-oo z(z - 1) 

The assertion follows now from Theorem A.l. □ 

Remark 4.1. 1. Using z ^-\) = zPT ~~ \ an< ^ substituting z — 1 for z we 
can write the variance-optimal initial capital for the European call option 
as 

Vo = 5o^w(log(|))-^(°)(log(| 
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with 

1 fR—j+ioo p zx 

:= / h(z + j) N — dz 

2m jR~j~ioo z 

in discrete time, respectively, 

1 rR-j+ico p rj(z+j)T+zx 

^0)( x ) := _/ dz 

2m jR-j-ioo z 

in continuous time. This resembles the pricing formulas for European calls 
in the Cox-Ross-Rubinstein, respectively, Black-Scholes, model. Note, how- 
ever, that tyki\x) may not be a distribution function in general. 

2. For the application of Lemma 4.1, we need slightly more than second 
moments of X\ and hence Sn- This seems counterintuitive because the payoff 
grows only linearly in Sn. It is, in fact, possible to derive the optimal hedge 
in the case where only second moments exist. The idea is to consider the 
difference of the call and the stock [cf. (33)]. Since the stock itself corresponds 
to the unit mass II = e\ , one immediately obtains an integral representation 
(2) of the call in the strip < K(z) < 1. In this case, II is a complex measure 
concentrated on {1} U (R + iR). 

Let us consider the Laplace transform representations of some simpler 
payoff functions. They are mostly taken from [42] and they can be derived 
by straightforward calculations from Theorem A.l. Interestingly, the put 
option payoff is expressed by the same integral as the call, but with the 
vertical line of integration to the left of zero, that is, 

i pR+ioo K^-~ z 

{K-s) + = — s z - -dz (R<0). 

2m jR-ioo z{z - 1) 

A related example is the payoff 

1 rR+ioo K^~ z 

(33) (s-K) + -s = — s z - -dz (0<R<1). 

2m JR-ioo z(z - 1) 

While this does not correspond to an option that arises in practice, it can be 
used to compute the variance-optimal hedge for calls and puts in a situation 
when the moment or cumulant function of the underlying asset exists in 
< 3i(z) < 2, but in no larger strip. This is actually the natural minimal 
integrability requirement in the present setup. 
The power call (cf. [43]) can be represented by 

1 rR+ico 2K 2 ~ Z 

«-*> + > = 2^ Loo 5 %(,-D(, -2) dZ <* >2 »' 
which generalizes to higher integer powers as 

i rR+ioo n^K n ~ z 

{{s-K)+r = — s* n , -dz (R>n) 

2m JR-ioo z(z - I) ■ ■ ■ (z - n) 
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and even to arbitrary powers a > 1 by 

1 rR+ioo 

((s - K) + ) a = / s z K a - z B(a + l,z-a)dz {R>a), 

2m jR-ioo 

where B denotes the Euler beta function, which can be expressed by the 
more familiar Euler gamma function 

The self-quanto call can be written as 

i rR+ioo K l ~ z 

(s-K) + s = / s z - -dz (R>2). 

K ' 2mJ R - ioo (z-l)(z-2) v ; 

The digital option with payoff function f(s) = Q{k,oo){ s ) coincides almost 
surely with the payoff function 

(34) f(s) = ^t {K} (s) + V,oc)(s) 

if the law of Sjv, respectively, St, has no atoms. Using statement 2 in The- 
orem A.l, the latter can be expressed as 

1 rR+ic zs-z 

f(s)= lim / s z dz (R>0). 

v ; c^oo 2vrz J R - ic z y ' 

This suggests applying the results of the previous sections to the measure 

1 K~ z 

(35) U(dz) = dz 

2m z 

in the case of the digital option. However, this measure is not of finite vari- 
ation. Nevertheless, the main statements remain valid if we interpret the 
integrals as Cauchy principal value integrals. 

Lemma 4.2. Theorems 2.1, 2.2 and 3.1-3.3 hold for the digital option 
(34) and the measure (35) if the integrals are interpreted in the principal 
value sense, that is, 

rR+ic 

H n :=P- lim S z h(z) N - n U(dz), 

c ^°° JR-ic 

rR+ic 



(36) £ n :=P-lim / S z n z\g{z)h{z) N ~ n Il{dz) , 

rR+ic rR+ic 

J := lim / / J (y,z)U(dy)U(dz) 
c ^°° JR-ic JR-ic 

and so forth, where P-lim refers to convergence in probability. In continuous 
time, the corresponding limit for £t(a>) is to be interpreted in (P(cLu) (g> dt) 
measure. 
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Proof. We will show the assertion in the continuous-time setting. The 
discrete-time case follows similarly. 

Step 1. For c e K+ define ffW := / C (5 T ) with 

pR+ic 

f(s):= s z U(dz). 

JR-ic 

Since K~ z jlixiz = — K~ z /2iriz, it follows that is real-valued. For s € 
M+, we have 



oo 2m \jR+ic \K J z JR-ic'\Kj z 

c'~*ooTTj c \ R + ix J 
The integrand equals 

s \ ( Rcos(xlog(s / K)) ^ R 2 sin(xlog(s / K)) sin(xlog(s/K)) 



K ) V R? + x 2 (R 2 + x 2 )x x 

Since sup cGR+ | / c °° sin j x ) dx\ < qq ( c f, [x], Section 5.2), it follows that 

sup \f(s)-f(s)\<us R 



for some u € K+. Consequently, (i?( c ) — H) 2 < u 2 S^ R £ L 1 for any c € M+, 
which implies that jj i n £ 2 by dominated convergence. 

Step 2. Denote by H = Hq + J Q T £ t <i5 t + Lt the Follmer-Schweizer decom- 
position of H, which exists, for example, by [38], Theorem 3.4. Moreover, 
let H[ c \ i[ c) and L ( t c) be defined as in Proposition 3.1 for the claim H^ c > . 
By Theorem 3.8 in [38], we have — > Ho, 

(37) E^j\^-^) 2 d{M,M) t ^0 

and E((L^ — Lt) 2 ) ^0 for c—> oo. Since and L are martingales, this 

implies — > in L 2 and hence in probability for any i £ [0, T] . Together 
with (14) we obtain 



f\d c) -Cu)dM u ^0, 
Jo 

[\^-L)dA u ^0 
Jo 



and hence 

rt _ 



Jo Jo 
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in probability for any t £ [0, T]. Moreover, we have C (c) £ 

m measure 

relative to P(du>) <g>dt [cf. (37) and (13)]. Together we obtain that Hq,^ and 
L 4 coincide with the expressions in Proposition 3.1 for H if the integrals are 
interpreted in the principal value sense. Theorems 3.1 and 3.3 now follow 
precisely as in Section 3. 

(c) ~ 

Step 3. Denote by Jq and Jo the variance of the hedging error for 
/j( c ) and H, respectively. In a Hilbert space, the mapping x i— > \\x — P(x)|| 2 
is continuous if P denotes the projection on some closed subspace. Hence 
Jq C ^ — > Jo for c — > oo. Since Theorem 3.2 is applicable to H^ c \ it follows that 
Jo coincides with Jo in (36). □ 

The log contract of [39] does not seem to fit into this framework because 
the logarithm has no Laplace transform. Nevertheless we can express it as a 
difference of two payoffs, namely its positive and negative part. The former 
has a Laplace transform for K(z) > and the latter for 3l(z) < 0, and we 
have 

1 rR+ioo ^ Y rR'+ioo ^ 

log(s) = 2^ L sZ ^ dz ~2ri L, s z # dz 

^I't-JR—ioc 6 Alii J R'— ioc £ 

with R' < and R > 0. In this case, II is a complex measure concentrated 
on (R' + iR)u(R + iM). 

Finally, let us emphasize again that the whole approach is linear in the 
claim. Hence, we immediately obtain the variance-optimal hedge for any 
linear combination of the payoffs above as, for example, bull and bear spreads 
and collars. 



5. Examples with numerical illustrations. In this section we illustrate 
how the approach is applied to concrete models that are considered in the 
literature. As an example, we provide numerical results for the normal in- 
verse Gaussian model. The other setups lead to similar figures. Recall that 
all quantities are discounted, as in the theoretical developments above. 

5.1. Discrete-time hedging in the Black-Scholes model. Suppose the un- 
derlying asset follows geometric Brownian motion with annual drift param- 
eter [i and volatility a. Then log returns per unit of time are normally 
distributed with mean // — <r 2 /2 and variance a 2 . 

Let us consider an option with maturity T and a positive integer N. If 
trading is restricted to times kT/N for k = 0, 1, . . . , N, the market becomes 
incomplete. Theorem 2.1 applies with the moment generating function 

m(z) = exp((( M -^), + ^)£). 
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If continuous trading is permitted, the Black-Scholes market is complete. 
Hence the hedging error is exactly zero. The variance-optimal capital and 
hedging strategy are given by the Black-Scholes price and delta hedging, 
respectively. It can be verified easily that this agrees in fact with the formulas 
in Theorems 3.1 and 3.2, where the relevant cumulant function is 



Note that the drift parameter [i affects the discrete-time results but has no 
influence in continuous time. 

5.2. Mertori 's jump diffusion with normal jumps. In the jump-diffusion 
model considered by [37], the logarithmic stock price is modeled as a Brown- 
ian motion with drift fi and volatility a plus occasional jumps from an inde- 
pendent compound Poisson process with intensity A. A particularly simple 
and popular case is obtained when the jumps are normally distributed, say 
with mean v and variance r. Then the moment function for discrete time is 



Note that Merton uses a slightly different parameterization. 

It turns out that the variance-optimal capital for a European call option 
may be negative for some parameters, for example, for continuous-time hedg- 
ing with So = 100, K = 110, T = 1 and fx = 0.01, a = 0.03, A = 0.01, v = 0.2, 
r = 0.02, we obtain Vq = —0.13. This illustrates, as mentioned above, that, 
in general, the variance-optimal initial capital is not a price. 

5.3. Hyperbolic, normal inverse Gaussian and variance gamma models. 
The hyperbolic, normal inverse Gaussian (NIG) and the variance gamma 
(VG) Levy processes are subfamilies or limiting cases of the class of general- 
ized hyperbolic models, which all fit in the general framework of this paper. 
We refer to [15] for further details. 

5.3.1. Hyperbolic model. The moment generating function in the hyper- 
bolic case is 



where K\ denotes the modified Bessel function of the third kind with index 1. 
Some care has to be taken if T/N is not an integer. The T/Nth power in 





and the cumulant function for continuous time is 
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(1) is, in fact, the distinguished T/Nth power (see [46], Section 7). The 
cumulant function equals 

KtjSy/aP -(J3 + zf) cllz \ 

Here Ln denotes the distinguished logarithm; see [46], Section 7. 

5.3.2. Normal inverse Gaussian model. The moment generating function 
of the normal inverse Gaussian model is given by 

m(z) = exp ((jjiz + 5(sja 2 -f3 2 - sj a 2 - ((3 + z) 2 ))^ 

Consequently, the cumulant function equals 

k(z) = fiz + 5(^Ja 2 -P 2 - ^Ja 2 - Q3 + z) 2 ). 

5.3.3. Variance gamma model. If we subordinate a Brownian motion 
with drift (5 by the gamma subordinator with parameters 5 and a, and 
add a linear drift with rate /x, we obtain a variance gamma Levy process. 
We refer to [35] for further details. 

The corresponding moment function for discrete time is 

m(z)=(e fiz ( " 0I ) S ) T/N . 

y ' \ \a-(3z-z 2 /2j ) 

The cumulant function needed for continuous-time hedging is given by 

k(z) =fiz + 5Ln( — — ). 

\a — pz — z A /2 J 

5.4. Numerical illustration. Figures 1-3 illustrate the results for a Euro- 
pean call in the normal inverse Gaussian model, compared to Black-Scholes 
as a benchmark. The time unit is a year. The parameters of the normal 
inverse Gaussian distribution are a = 75.49, [3 = —4.089, 5 = 3.024 and 
(i = —0.04, which correspond to annualized values of the daily estimates 
from [45] for Deutsche Bank, assuming 252 trading days per year and dis- 
counting by an annual riskless interest rate of 4%. 

The parameters for the benchmark Gaussian model are chosen such that 
both models lead to returns of the same mean and variance. We consider a 
European call option with strike price K = 99 and maturity T = 0.25, that 
is, three months from now. 

Figure 1 shows the variance-optimal initial capital as a function of the 
stock price in the NIG model for both continuous and discrete time with 
N = 12 trading dates, that is, weekly rebalancing of the hedging portfolio. 



k(z) = Ln 
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The Black-Scholes price is plotted as well for comparison. One may observe 
that the three curves cannot be distinguished by eye, that is, they do not 
differ much in absolute terms. A similar picture is obtained in Figure 2 for 
the initial hedge ratio as a function of the initial stock price. The Black- 
Scholes delta provides a good proxy for the optimal hedge in the NIG model 
for both continuous and weekly rebalancing. As a result one may say that 
the Black-Scholes approach produces a reasonable hedge for the European 
call even if real data follow this rather different jump-type model. 

The similarity ceases to hold when it comes to the hedging error, which 
vanishes in a true Black-Scholes world. Figure 3 shows the variance of the 
hedging error as a function of the number of trades, going from N = 1 (static 
hedging) to N = 63 (daily rebalancing). The horizontal line in Figure 3 
indicates the variance of the hedging error for continuous rebalancing in the 
NIG model. The two decreasing curves refer to the discrete hedging error in 
the NIG and the Gaussian case, respectively. In the latter case it converges 
to 0, which is the error in the limiting Black-Scholes model. As far as the 
size is concerned, the variance of the error for N = 12 trading dates (weekly 
rebalancing) in the NIG setup (1.04 ~ 1.02 2 ) equals roughly the sum of the 
error in the corresponding Gaussian model (0.83, due to discrete rather than 
perfect hedging) and the inherent error in the continuous-time NIG model 
(0.257, due to incompleteness from jumps). The standard deviation 1.02 of 

25 1 1 1 1 1 1 1 1 



Discrete NIG endowment 
Continuous NIG endowment 
Black-Scholes price 
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Fig. 1. Variance- optimal initial capital for normal inverse Gaussian returns. 




Fig. 2. Variance- optimal initial hedge for normal inverse Gaussian returns. 
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Fig. 3. Variance of the hedging error for normal inverse Gaussian returns. 
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the hedging error in the discrete NIG case may be compared to the Black- 
Scholes price 4.50 of the option. 

APPENDIX: BILATERAL LAPLACE TRANSFORMS 

Definition A.l. Let f:R^C be a measurable function. The (bilat- 
eral) Laplace transform f is given by 

(A.l) /(*)=/ f(x)e~ zx dx 

J —oo 

for any z € C such that the integral exists. 

In the literature, the Laplace transform / is also denoted by C[f(x);z] or 
by Cn[f(x);z] when it is necessary to distinguish the bilateral from the usual 
(unilateral) Laplace transform. The latter is defined by the same integral, 
but starting from instead of — oo. 

We say that the Laplace transform f(z) exists if the Laplace transform 
integral (A.l) converges absolutely, or, in other words, if it exists as a proper 
Lebesgue integral as opposed to an improper integral. The following lemma 
shows that the domain of a Laplace transform is always a vertical strip in 
the complex plane. It may be empty, degenerate to a vertical line, a closed 
or open left or right half-plane, or all of C. 

Lemma A.l. Suppose that f(a) and f(b) exist for real numbers a <b. 
Then f(z) exists for any z € C with a < ^t(z) < b. 

Proof. The proof is obvious because \ f(x)e~ zx \ = \f(x)\e~^ z " >x < \f(x) x 
e- ax \ + \f(x)e- bx \. □ 

From 

„ r+oo r+oo 

(A.2) f(u + iv)= / f(x)e- {u+lv)x dx= / e ux f(-x)e lxv dx, 

J —oo J —oo 

we see that C[f(x); u + iv] = T\e ux f{— x)\ v], where the last expression de- 
notes the Fourier transform of the function x \— > e ux f(—x). Hence all prop- 
erties of the bilateral Laplace transform can be reformulated in terms of the 
Fourier transform and vice versa. 

There are many inversion formulas for the Laplace transform known in 
the literature. We will use the so-called Bromwich inversion integral, which 
can be justified by the following theorem. 

Theorem A.l. Suppose that the Laplace transform f(R) exists for R £ 

E. 
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1. If v i— > f(R + iv) is integrable, then x f{x) is continuous and 

1 rR+ioo _ 

f(x) = — f(z)e zx dz forxtR. 

2m jR-ioo 

2. If f is of finite variation on any compact interval, then 

1 1 rR+ic „ 

Iim-(/(x + e) + /(a;-e))= lim / f(z)e zx dz forx£R. 

e^O 2 c->oo 2m JR-ic 

Proof. The first statement follows from [44], Theorem 9.11 and (A. 2). 
For the second assertion, see [10], Satz 4.4.1. □ 
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This paper is an extended and revised version of [27]. 
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